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Abstract 

The sidereal time dependence of MiniBooNE v e and v e appearance data are 
analyzed to search for evidence of Lorentz and CPT violation. An unbinned 
Kolmogorov-Smirnov test shows both the v e and P e appearance data are com- 
patible with the null sidereal variation hypothesis to more than 5%. Using 
an unbinned likelihood fit with a Lorentz- violating oscillation model derived 
from the Standard Model Extension (SME) to describe any excess events 
over background, we find that the v e appearance data prefer a sidereal time- 
independent solution, and the v e appearance data slightly prefer a sidereal 
time-dependent solution. Limits of order 10~ 20 GeV are placed on combina- 
tions of SME coefficients. These limits give the best limits on certain SME 
coefficients for — > v e and — > v e oscillations. The fit values and limits of 
combinations of SME coefficients are provided. 

Keywords: MiniBooNE Neutrino oscillation Lorentz violation PACS: 
11.30.Cp 14.60.Pq 14.60. St 



1. Introduction of Loremtz violation 

Violation of Lorentz invariance and CPT symmetry is a predicted phe- 
nomenon of Planck scale physics, especially with a spontaneous violation [1], 
and it does not require any modifications in quantum field theory or general 
relativity. Since neutrino oscillation experiments are natural interferome- 
ters, they can serve as sensitive probes of space-time structure. Neutrino 
oscillations have the potential to provide the first experimental evidence for 
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Lorentz and CPT violation through evidence of oscillations that deviate from 
the standard L/E dependence [2], or that show sidereal time dependent os- 
cillations as a consequence of a preferred direction in the universe [3]. 

In this paper, we test the MiniBooNE — > v e and — > v e oscillation 
data [4, 5] for the presence of a Lorentz violation signal. Similar analyses 
have been performed in other oscillation experiments, including LSND [6], 
MINOS [7], and IceCube [8]. Naively, experiments with longer baselines 
and higher energy neutrinos would be expected to have better sensitivity to 
Lorentz violation, because small Lorentz violating terms are more prominent 
at high energy, where neutrino mass terms are negligible. However, some 
Lorentz violating neutrino oscillation models mimic the standard massive 
neutrino oscillation energy dependence [9]. Then, in this case, the signal 
may only be seen in sidereal variations of oscillation experiments. 

2. MiniBooNE experiment 

MiniBooNE is a v e (u e ) appearance short baseline neutrino oscillation 
experiment at Fermilab. Neutrinos are created by the Booster Neutrino 
Beamline (BNB), which produces a 93% (83%) pure (z/ M ) beam in neutrino 
(anti-neutrino) mode, determined by the polarity of the magnetic focusing 
horn. The MiniBooNE Cherenkov detector, a 12.2 m diameter sphere filled 
with mineral oil, is used to detect charged particles from neutrino interactions 
and is located 541 m from the neutrino production target. It is equipped with 
1,280 8" PMTs in an optically separated inner volume, and 240 8" veto PMTs 
in an outer veto region. Details of the detector and the BNB can be found 
elsewhere [10, 11]. Charged leptons created by neutrino interactions in the 
detector produce Cherenkov photons, which are used to reconstruct charged 
particle tracks [12]. The measured angle and kinetic energy of charged leptons 
from neutrino interactions are used to reconstruct the neutrino energy, E® E , 
for each event, under the assumption that the target nucleon is at rest inside 
the nucleus and the interaction type is a charged current quasielastic (CCQE) 
interaction [13]. 

For this analysis, we use the background and error estimates from [14] 
(neutrino mode) and [15] (antineutrino mode). For neutrino mode, data from 
6.46 x 10 20 protons on target (POT) are used. An excess in the "low energy" 
region (200< i^- E (MeV)<475) was observed, with 544 events reported as 
compared to the prediction, 409.8 ± 23.3(stat.)±38.3(syst.). Interestingly, 
this excess does not show the expected L/E energy dependence of a simple 
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two massive neutrino oscillation model. Additionally, it is not consistent with 
the energy region expected for the "LSND" signal [16]. For the anti-neutrino 
mode analysis (5.66 x 10 20 POT), MiniBooNE observed a small excess in the 
low energy region, and an excess in the region 475< -E^ £ (MeV)<1300. The 
excess in this "high energy" region is found to be consistent with the LSND 
signal, assuming a two massive neutrino hypothesis, but remains statistically 
marginal. In the "combined" region (200< £^ E (MeV)<1300), MiniBooNE 
observed 241 V e candidate events as compared to the prediction, 200.7 ± 
15.5(stat.)±14.3(syst.). 

Although the conflict between MiniBooNE neutrino and anti-neutrino 
mode results can be resolved in models without CPT violation [17], CPT vi- 
olation is a viable option. Since CPT violation necessarily implies violation 
of Lorentz invariance within interactive quantum field theory [18], we are 
in a well-motivated position to search for Lorentz and CPT violation using 
the MiniBooNE data. In fact, proposed models motivated by Lorentz viola- 
tion [19, 20] can already accommodate world data including the MiniBooNE 
and LSND excesses with a small number of free parameters. Evidence for 
sidereal variation in the MiniBooNE excesses would provide a distinctive 
direct signal of Lorentz violation. 

3. Analysis 

We use the Standard Model Extension (SME) formalism for the general 
search for Lorentz violation [21]. The SME is an effective quantum field 
theory, and the minimum extension of the Standard Model including particle 
Lorentz and CPT violation [21]. A variety of data have been analyzed under 
this formalism [22], including neutrino oscillations [6, 7, 8]. In the SME 
formalism for neutrinos, the evolution of a neutrino can be described by an 
effective Hamiltonian [3], 



Here, E and are the energy and the 4-momentum of a neutrino, and (ol)^, 
and (cl) 1 ^ are CPT-odd and CPT-even SME coefficients in the flavor basis. 
Under the assumption that the baseline is short compared to the oscillation 
length [23], the — > u e oscillation probability takes the form, 




(1) 
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(he) 2 



|(C) eM + (A s ) en sinw ffi T ffi + (A c )en coscj e T e 
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+(£ s ) eM sin2w ffi Te + {B c ) efl cos 2w ffi T ffi | 2 . (2) 

This probability is a function of sidereal time T e . Four parameters, (A s ) efl , 
(Ac) ef i, (Bs)e(i, an d (B c )en are sidereal time dependent, and (C) eil is a sidereal 
time independent parameter. We use a baseline distance of L =522.6 m, 
where the average pion decay length is subtracted from the distance between 
the neutrino production target and detector. And u ffi is the sidereal time 
angular frequency described shortly. These parameters are expressed in terms 
of SME coefficients and directional factors [23]. The same formula describes 
the Up —¥ v e oscillation probability by switching the signs of the CPT-odd 
SME coefficients. We neglect the standard neutrino mass term, to 2 / 'E <C 
10 -20 GeV, which is well below our sensitivity which is discussed later. 

For this analysis, we convert the standard GPS time stamp for each event 
to local solar time (period 86400.0 sec) and sidereal time (period 86164.1 sec). 
We then define the local solar time angular frequency u Q = 86 2 q Q Q (rad/s) 
and the sidereal time angular frequency w e = 86 2 g 4 1 (rad/s). The time origin 
can be arbitrary, but we follow the standard convention with a Sun-centered 
coordinate system [6]. We choose a time- zero of 58 min after an autumnal 
equinox of 2002 (September 23, 04:55GMT), so that this serves not only 
as the sidereal time-zero, but also as the solar time-zero since Fermilab is 
on the midnight point at this time. The local coordinates of the BNB are 
specified by three angles [23], colatitude x — 48.2°, polar angle 6 = 89.8°, 
and azimuthal angle <fi = 180.0°. 

Any time dependent background variation, such as the time variation 
of detector and BNB systematics, are important. To evaluate these, we 
use our high statistics CCQE samples (Figure 1). These data are from our 
z/^CCQE double differential cross section measurement sample [24] composed 
of 146, 070 events (5.58 x 10 20 POT), and our P^CCQE candidate sample [25] 
composed of 47,466 events (5.66 x 10 20 POT). The ^(P^CCQE local solar 
time distribution exhibits ±6(3)% variation. 

The same variation in local solar time is observed in the POT for both 
neutrino and anti-neutrino mode data taking periods. Therefore, the POT 
variation is the dominant time dependent systematic error. The amplitude is 
negligible in z/^CCQE sidereal time distribution, however, it persists in ~3% 
variations in z/^CCQE sidereal time distribution. We evaluate the impact of 
this variation on our analysis by correcting POT variation event by event 
in v e (p e ) candidate data. It turned out the correction only has a negligible 
effect. Thus we decided to use unweighted events. This also simplifies the 
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Figure 1: The top (bottom) plot shows the f^CCQE (P^CCQE) local solar time distri- 
bution. The solid curves are fit functions extracted from the CCQE event distributions, 
and the dotted curves are from the POT distributions (arbitrary units) during the same 
period of data taking. The dashed line shows a flat distribution. 

unbinned likelihood function used in later analysis. Figures 2 and 3 show 
the v e and v e oscillation candidates sidereal time distributions both with 
and without the POT correction. These plots verify that time- dependent 
systematics are negligible in this analysis. 

To check for a general deviation from a flat distribution (null sidereal vari- 
ation hypothesis), we perform an unbinned Kolmogorov-Smirnov test (K-S 
test) [26] as a statistical null hypothesis test for both the v e and v e sam- 
ples. The K-S test is suitable in our case because it is sensitive to runs in 
distributions, which may be a characteristic feature of the sidereal time de- 
pendent hypothesis. Table 1 gives the result. The K-S test is applied to 
the low energy, high energy and combined regions, for both neutrino, and 
anti-neutrino mode data. To investigate the time dependent systematics, we 
also apply the K-S test to the local solar time distribution. The test shows 
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low energy high energy 


combined 




solar sidereal solar sidereal 


solar sidereal 




Neutrino mode 




<E V > 


0.36 GeV 0.82 GeV 


0.71 GeV 


#evt 


544 420 


964 


P(KS) 


0.42 0.13 0.81 0.64 


0.64 0.14 




Anti-neutrino mode 




<E V > 


0.34 GeV 0.78 GeV 


0.60 GeV 


#evt 


119 122 


241 


P(KS) 


0.62 0.15 0.79 0.39 


0.69 0.08 



Table 1: A summary of K-S test results on the sidereal and local solar time distributions. 
The top table is for v e candidate data, and the bottom table is for v e candidate data. The 
three rows show the average neutrino energy of each sample, number of events, and the 
K-S test compatibility with the null hypothesis. The test is performed in three energy 
regions, and for both solar local time and sidereal time distributions. 

none of the twelve samples has less than 5% compatibility (~ 2a), which 
we chose as a benchmark prior to the analysis. Hence, all samples are com- 
patible with the null sidereal variation hypothesis. Interestingly, the sidereal 
time distributions tend to show lower compatibility with a flat hypothesis, 
but not by a statistically significant amount. These results indicate that any 
sidereal variation extracted from our data, discussed below, is not expected 
to be statistically significant. 

To fit the data with the sidereal time-dependent model, we use a gen- 
eralized unbinned maximum likelihood method [27]. This method finds the 
best fit model parameters by fitting data with a log likelihood function I. It 
is suitable for our analysis because this method has the highest statistical 
power for a low statistics sample. In this method, the log likelihood function 
i is constructed by adding 1^ from each event. After dropping all constants, 
£i has the following expression, 

Here, A" is the number of observed candidate events, fi s is the predicted num- 
ber of signal events which is given by the time integral of Eq. 2 together with 
the estimated efficiency, \x\, is the predicted number of background events, 
T s is the probability density function (PDF) for the signal and is a function 
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of sidereal time and the fitting parameters (Eq. 2 with proper normaliza- 
tion), T\y is the PDF for the background, is the la error on the predicted 
background, and m is the central value of the predicted total background 
events. Two sources contribute equally to the background: intrinsic beam 
background and mis-identification (mainly ir°s). Their total variation is as- 
signed as the systematic error, assumed to be time independent. Details can 
be found in [4, 5]. The parameter space is scanned (grid search method) 
to find the largest £, or the maximum log likelihood (MLL) point, and this 
MLL point provides the combination of the best fit (BF) parameters. The log 
likelihood function includes six parameters, five that are functions of SME 
coefficients, and one for the background. However the background term is 
constrained within a ±la range. Neither the neutrino nor the anti-neutrino 
mode data allow us to extract errors if we fit all five parameters at once, due 
to the high correlation of parameters. Therefore, we set (B s ) efl and (B c ) efl to 
zero and concentrate on three parameters ((C) e ^, (A s )en, and (A c ) eiJl ) for the 
uncertainty estimates. Since the five parameter fit is quantitatively similar 
to the three parameter fit, we will focus the discussion on the results on the 
three parameter fits. This three parameter fit also corresponds to the case 
with only CPT-odd SME coefficients [23]. 

4. Results 

Figure 2 shows the neutrino mode low energy region fit results. The top 
three plots show the three projections of three dimensional parameter space. 
Because of the square of fitting parameters in the PDF, the BF point has a 
sign ambiguity and is always duplicated. The la and 2a contours are formed 
from a constant slice of the log likelihood function in the three dimensional 
parameter space. To avoid under coverage, these slices are expanded until 
they enclose 68% (la) or 95% (2a) of BF points for the three parameter fit 
of simulated, or "fake" , data sets with the signal. Note that because fitting 
parameters are not linear in the PDF, twice the la error does not yield the 
2a error. 

A null sidereal variation hypothesis, or a flat solution, is equivalent to a 
three or five parameter fit solution where only the (C) eM parameter is nonzero. 
The fit to neutrino data favors a nonzero solution only for the (C) efM term. 
The bottom plot in Figure 2 shows data plotted against curves corresponding 
to the flat solution and the best fits for three and five parameter functions. 
Since all three curves are close to each other, the solution of neutrino mode is 
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v— mode BF 2a limit v— mode BF 



2cr limit 



|(C) e(U | 3.1 ±0.6 ±0.9 <4.2 0.1 ±0.8 ±0.1 < 2.6 
\(As)eiA 0.6 ±0.9 ±0.3 <3.3 2.4 ±1.3 ±0.5 < 3.9 
\{Ac)e^\ 0.4 ±0.9 ±0.4 <4.0 2.1 ±1.2 ±0.4 < 3.7 

SME coefficients combination (unit 10~ 2 ° GeV) 
WU\ ±[(«l)^ +0.75(a L )fJ- < E > [1.22( Ci )^' + 1.50(c L )^ + 0.34(c L )^] 

±[0.66(a i )D- < £ > [1.33( Ci )ST + 0.99(c L )r/] 
\(AcU\ ±[0.66(a L )*]- < E > [1.33(c L )^f + 0.99(c L )£ z ] 



Table 2: The fit parameters for the neutrino mode low energy region and the anti- neutrino 
mode combined region. The BF points are the MLL points of the log likelihood function, 
here top rows from left to right, BF values, la statistical, and systematic errors. The 
2a limits are also shown. Bottom rows show detailed expressions of each sidereal fit 
parameter in terms of SME parameters, and directional factors [23]. The upper (lower) 
sign of (at)^ terms is applied for neutrino (anti-neutrino) results, due to the CPT-odd 
nature. The average neutrino energy "< E >" is 0.36 GeV for the neutrino mode low 
energy region and 0.60 GeV for the anti-neutrino mode combined region (Tab. 1). 

dominated by the sidereal time- independent component. To find the signifi- 
cance of time dependence over the flat distribution, fake data sets without a 
signal are formed where the v e candidate events are simulated without any 
time structure. The MLL difference between the three parameter fit and 
the flat solution is used to form a A% 2 , and the expected Ax 2 distribution 
is determined by testing 500 random distributions from the fake data sets. 
This test shows that there is a 26.9% chance that a random distribution of 
v e candidate events would yield a A% 2 value equal to, or greater than, the 
value observed for the data. This result is consistent with the sensitivity of 
this experiment. We estimate our sensitivity to the limitted case. First, a 
2a threshold is set from this A% 2 distribution. Then, time dependent ampli- 
tudes were incrementally increased in the model until the 2a threshold was 
exceeded. When we assume (C) efM ^ and (A s ) efl = (A c ) efM ^ the 2a dis- 
covery threshold of sidereal time dependent amplitudes from v e {p e ) candidate 
data statistics are 0.8(1.1) xl0~ 20 GeV. 

Figure 3 shows the analogous fit results for the anti-neutrino mode com- 
bined energy region. Due to lower statistics, the combined region is used 
rather than dividing the data into two subsets. Unlike the neutrino mode 
low-energy region, the (C) eAt parameter no longer significantly deviates from 
zero. The fit to anti-neutrino data favors a non-zero solution for the (A s ) e ^ 
and (Ac)e[i parameters at the nearly 2a level. Performing the same A% 2 test 
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Figure 2: (color online) Three parameter fit results for the neutrino mode low energy 
region. The top three plots show the projection of three dimensional parameter space. 
The dark (light) shaded area shows the la (2a) contour in each projection. The stars 
show the BF points. The bottom plot shows the curves corresponding to the flat solution 
(dotted), three parameter fit (solid), and five parameter fit (dash-dotted), together with 
binned data (solid marker). The POT corrected data are also shown in open circle marker. 
Here, the fitted background is shown as a dashed line, and the BF value is 1.00 (i.e., 
equivalent to the central value of the predicted background). 

as is outlined above results in only 3.0% of the random distributions from 
the v e candidate events having a A% 2 value exceeding the value observed for 
the data. Note that this is consistent with the 8% compatibility with a flat 
hypothesis found with the K-S test (Tab. 1). 

Table 2 shows fit parameters for the neutrino mode low energy region 
and anti-neutrino mode combined region. All errors are estimated from la 
contours of parameter space projections. Errors are asymmetric, but we 
choose the larger excursions as the symmetric errors for simplicity. The 2a 
contours provide the limits. In principle, these fit parameters are complex 
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Figure 3: Three parameter fit results for the anti-neutrino mode combined region. Nota- 
tions are the same as Fig. 2. Here, the BF value for the fitted background is 0.97 (3% 
lower than the central value of the predicted background) . 

numbers. Here, all parameters are assumed to be real. A naive estimation 
from Tab. 2 indicates possible SME coefficients to satisfy the MiniBooNE 
data are of order 1CT 20 GeV (CPT-odd), and lCT 20 to 1CT 19 (CPT-even). 
However, these SME coefficients are too small to produce a visible effect for 
LSND [6]. On the other hand, any SME coefficients extracted from LSND [6] 
predict too large of a signal for MiniBooNE. Therefore, a simple picture 
using Lorentz violation to explain both data sets leaves some tension, and a 
mechanism to cancel the Lorentz violating effect at high energy [3, 19, 20] is 
needed. 

5. summary 

In summary, we performed a sidereal time variation analysis for Mini- 
BooNE v e and v e appearance candidate data. For the neutrino mode low 
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energy region, K-S test statistics indicate the null hypothesis is compatible 
at the 13% level, and the relative improvement in the likelihood between 
the null hypothesis and the three parameter fit occurs 26.9% of the time in 
random distributions from a null hypothesis. Analysis of the combined en- 
ergy region in anti-neutrino mode results in a K-S test that indicates a 8% 
compatibility with the null hypothesis, however the relative improvement in 
the likelihood between the null hypothesis and the three parameter fit only 
occurs 3.0% of the time in random distributions from a null hypothesis. The 
limits of fit parameters, 1CT 20 GeV, are consistent with Planck scale sup- 
pressed physics. This is the first sidereal variation test for an anti-neutrino 
beam of ~1 GeV energy and ~500 m base line. These limits are currently 
the best limits on the sidereal-time independent (a L ) e ^ and (cl) 6M SME co- 
efficients. These limits can be significantly improved by long baseline v e (P e ) 
appearance experiments, such as T2K [28] and NOvA [29]. 
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